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We report on the measurement of the neutron radiation hardness of silicon photomultipliers (SiPMs) manufactured by Hamamatsu 
Corporation in Japan and SensL in Ireland. Samples from both companies were irradiated by neutrons created by a 1 GeV electron 
beam hitting a thin lead target at Jefferson Lab Hall A. More tests regarding the temperature dependence of the neutron radiation 
damage and self-annealing were performed on Hamamatsu SiPMs using a calibrated Am-Be neutron source from the Jefferson Lab 
Radiation Control group. As the result of irradiation both dark current and dark rate increase linearly as a function of the 1 MeV 
equivalent neutron fluence and a temperature dependent self-annealing effect is observed. 
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1. Introduction : 

A Silicon photomultiplier (SiPM) is a photon-counting de- . 
vice consisting of multiple avalanche photodiode (APD) pixels . 
operating in Geiger mode. It is also known as Multi-Pixel Pho- . 
ton Counter (MPPC). Each APD pixel of the SiPM outputs a 
pulse signal when it detects photons and the signal output from 
the SiPM is the total sum of the signals from all APD pixels. 
The SiPM offers the high performance needed in photon count- 
ing and is used in diverse applications for detecting extremely 
weak light intensities at the photon-counting level. 

For Hall D At Jefferson Lab [ 1|, we investigated the use of 
SiPMs to collect light from the Barrel Calorimeter. One of the 
requirements for such devices is sufficient radiation hardness to 
withstand many years of operation. As the neutron background 
is expected to be the major source of radiation damage in the 
HallQ we did a series of tests of the neutron radiation damage to 
SiPMs at various conditions to evaluate the life time of SiPMs 
in Hall D. 



2. Radiation Damage in Silicon Detectors 

The bulk damage in silicon detectors caused by hadrons or 
higher energy leptons and photons is primarily due to displace- 
ment of primary knock-on atoms from the lattice |j2)- For 
neutrons or electrons with kinetic energy above 175 eV and 29 
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260 keV, respectively, they will start to generate Frenkel pairs 
(a pair of a silicon interstitial and a vacancy) along their tracks 
in silicon material. With higher energy, more than 35 keV for 
neutrons and 8 MeV for electrons, a dense cluster of defects 
will be formed at the end of the primary PKA track. 
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Figure 1 : Effective damage to Silicon detectors relative to 1 MeV neutron. Data 
at different energy ranges are taken from References 1 3 4 5 |. 

The defects will effect silicon detector's performance in var- 
ious levels depending on their concentration, energy level and 
the respective electron and hole capture cross-section. For in- 
stance, defects with energy levels in the middle of the forbidden 
gap acting as recombination/generation centers are responsible 
for an increase of the reverse current. Interactions with dopants 
change the effective doping concentration and therefore change 
the operating voltage of the detector. Finally, defects acting as 
trapping centers reduce the charge collection efficiency. 

The radiation damage of neutrons to Silicon detectors has 
been extensively summarized in the literature. In this paper, we 
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took data from References J3j |4j [5l to convert the neutron flux 
to the 1 MeV neutron fluence based on the effective damage 
caused by neutrons with different energies. The weight of such 
a conversion is plotted in Fig.[T] 

3. Test with Electron Beam 

In addition to the goal of testing the neutron radiation dam- 
age to SiPMs, we also confirmed our knowledge of the neutron 
background. 



3.1. Test Setup 




Figure 2: A 4 X 4 SiPM array of 3 X 3 mm SiPM cells from Hamamatsu and 
individual cells are marked by dashed lines. 

The irradiation was carried out during the PRex experi- 
ment [6 1 at Jefferson Lab Hall A, with a 1 GeV electron beam 
incident on a 0.5 mm Pb target. During the two-day test, one 
Hamamatsu unr0and one SensL unij^Jwere irradiated. Both 
units are 4x4 arrays of 3 x 3 mm 2 SiPMs and were powered 
to a gain of 0.75xl0 6 . A photograph of the Hamamatsu SiPM 
array is shown in Fig. [2] 

Both SiPM units were placed inside a dark box together with 
their pre-amplifiers. The light from a pulsed LED was guided 
into the box through an optical fiber and then diffused by a 
diffuser to provide uniform light on both units. The box was 
positioned 20 meters away from the Hall A Pb target and 135 
degrees backwards to the beam direction as illustrated in Fig. [3] 
to reduce the effect from other sources of radiation such as pho- 
tons and charged particles. While the box had direct view of 
the target, the rest of the equipment was shielded by a con- 
crete wall. The shape of the output signals, including amplitude 
and width, was continuously recorded by an oscilloscope. In 



2 A preproduction unit of SI 0943-025 8(X) MPPC with 50 /im pixels, equiv- 
alent to new SI 2045. 

3 A SPMArray unit based on ceramic design with 35 lira pixels - Wafer Batch 76 
Code X4151-05 using SPM3035 design. 77 
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Figure 3: The floor plan and readout scheme for the SiPM radiation test in 
Hall A. 

the middle of the irradiation period, the power supplies of both 
SiPMs were turned off intentionally to see whether the power- 
ing condition affects the radiation damage. After the irradiation, 
both SiPMs were stored at room temperature, 20-25°C, for the 
self-annealing test. 
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Figure 4: The simulated neutron energy spectrum at the SiPM testing area in 
Hall A with 50 /jA 1 GeV electron on a 0.5 mm Pb target. 

A BF3 neutron probe was positioned next to the dark box to 
monitor the relative change of the neutron flux in real time. In 
order to obtain the absolute reading of the effective 1 MeV neu- 
tron fluence, a couple of the same type of Hamamatsu SiPMs 
were later irradiated by a calibrated AmBe neutron source to a 
similar damage level. As the AmBe source has a well known 



narrow energy spectrum peaking at about 4 MeV (T), its flu- 
ence was calculated by convoluting the neutron flux spectrum 
with the effective damage curve shown in Fig.[T] Then the flu- 
ence in Hall A was calculated by comparing the damage the 
SiPMs received in both cases. 

As a result, we determined that the two SiPMs in Hall A re- 
ceived a fluence of about 3.7 x 10 9 n eq /cm 2 . The fluence mea- 
surement also provides a good bench mark of our knowledge of 
the radiation level in the experimental halls. The neutron flux in 
Hall A was simulated in a GEANT3 framework customized to 
the electron-beam environment [8]|9][10). The resulting energy 
spectrum is shown in Fig.|4]and the fluence obtained from such 
a simulation is consistent with the measured value within 50%. 
The same code predicts that the 3.7 x 10 9 n eq /cm 2 fluence will 
be reached in about 13 years in Hall D with its high intensity 
GlueX running^Jon a 30 cm liquid Hydrogen target. 
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94 3.2. Results 



Figure 6: (Color online) Signals from the Hamamatsu SiPM array recorded by 
the oscilloscope during the irradiation. 



95 3.2.1. Dark Current 
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Figure 5 : (Color online) The increase of the dark current of SiPMs as a function 
of the neutron fluence during the two-day irradiation in Hall A. The curves are 
fits using second order polynomials. The period with no data marked by grey 
dash lines corresponds to the time when the sensors were not powered, but 
continued to be irradiated. 



The change of the SiPM dark current as a function of the ac- 
cumulated neutron fluence is plotted in Fig. [5] As the beam was 
turned on, the dark current of both SiPMs started to increase im- 
mediately. By comparing the trends of the damage before and 
after the period when the power was turned off, one can see 
that the neutron damage remains the same no matter whether 
the unit is powered or not. Over the course of the test, the dark 
current increased by a factor of about 10 for the SensL SiPM 
Array, 160 /vA -> 1.6 mA, and 25 for the Hamamatsu SiPM 
An-ay, 8 fiA -> 200 fiA. 
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Figure 7: (Color online) The impact of neutron radiation damage on the sig- 
nal shape, amplitude and width, of the SiPM output. The amplitude dropped 
slightly while no change was observed in the width. The curves are fits using 
second order polynomials. 
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3.2.2. Signal 

The output signals from both SiPMs stayed relatively sta- 
ble in contrast to the dramatic change of the dark current. The 
amplitude and width (50% to 50%) are plotted in Fig. [6] and 
[7] While the width shows no noticeable change, the amplitude 
dropped by about 10%. 

3.2.3. I-V Curve 
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Figure 8: (Color online) The I-V curve of Hamamatsu SiPM array before and 138 
after the irradiation. The dark current after irradiation is scaled by a 0.04 for a 
better visual comparison and clearly the break down voltage was not effected 139 
by the neutron radiation. 14 " 

141 

113 The current vs. voltage (TV) curves before and after irradi- 142 

114 ation are also compared for both units, and the comparison of 

115 the Hamamatsu SiPM array is shown in Fig. [8] Other than ani43 
us overall change in scale, the I-V curves stay the same for both 

144 

117 units and indicate that the break down voltage of the SiPM is 

° 145 

118 not impacted by the neutron radiation. 

146 

119 3.2.4. Self-Annealing 

148 

120 Following the delivery of an intense prompt dose of radia- 

121 tion, one finds that not all damage to the lattice is permanent. 

122 So we studied the response of the sensors as they rested with- 

151 

123 out further irradiation (self-annealing). The results of the self- 

124 annealing at room temperature after the prompt irradiation are 

125 plotted in Fig. [9] and 10 Both the dark current and signal am- 

126 plitude recovered over time with a time constant close to 10 

127 days. About half of the damage to the dark current recovered 

128 and the signal amplitude completely returned to the level before 

129 the irradiation. 

158 
159 

130 4. Temperature Tests with a AmBe Neutron Source , 60 

131 The Hall A irradiation test revealed that the lifetime of the 61 

132 SiPM will be marginal in Hall D given the experimental re- 

133 quire ment on the dark rate. Fortunately, a lower dark rate at 163 
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Such a high intensity refers to a Bremsstrahlung photon flux of about 11 
1 00 MHz/GeV close to the 1 2 GeV endpoint. 1 1 



Figure 9: (Color online) The decline of dark current during SiPM's self- 
annealing at room temperature after the irradiation. The time constant for both 
samples is about 10 days and approximately half of the damage recovered. The 
data are fitted by an exponential function described in Eq. 



a fixed gain can be achieved if the SiPM is cooled. Cooling 
the SiPMs to 5°C will reduce the dark rate to about 1/3 com- 
pared with 20°C and will certainly allow more room for the 
increase of the dark rate caused by the neutron radiation dam- 
age. However, whether such a dependence will be effected by 
radiation damage was unknown, therefore a systematic study 
of the temperature dependence of the neutron radiation damage 
and annealing was performed using a calibrated AmBe neutron 
source. 

4.1. Test Procedure 

Twelve lxl mm 2 SiPM units from Hamamatsu with 50 /vm 
pixels (PN# S10362-11-050C) were irradiated during this test. 

Since the previous Hall A irradiation test shows no correla- 
tion between the damage and the powering condition, all the 
SiPMs were not powered during the irradiation or annealing 
except when their dark currents were measured. If not specifi- 
cally mentioned, the dark current was always measured at room 
temperature regardless of the SiPM's irradiation or annealing 
temperatures. This allows direct comparison of results between 
samples regardless of the temperature at which they were irra- 
diated or annealed. It is assumed that taking the test samples 
to room temperature during the short time of the measurement 
does not significantly influence the results. The unit was then 
powered off and put back to its previous temperature after the 
dark current was measured. All the SiPMs were powered to a 
gain of 0.84 x 10 6 during the dark current measurements. The 
gain was determined using the ADC spectra and more details 
can be found in Sec. 14.2.21 

The test consists of the following steps as shown in Fig. [TT[ 

1. In the first stage, six units were irradiated at -5°C while 
the remaining six units were irradiated at room temperature, 
~ 25°C. The irradiation lasted four days and the total fluence 
each unit received was about 1.4 x 10 9 n eq /cm 2 . 
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Such a 110 nA difference suggests two possibilities, one is a 
temperature dependence of the radiation damage and the other 
one is a temperature dependence of the damage recovery. 

For the units annealed at 60° C during the first annealing, the 
average dark currents of the units from the -5°C and 25°C ir- 
radiation groups dropped to 365 + 61 nA and 341 ±23 nA, 
respectively. The consistency of these two values excludes the 
temperature dependence of the radiation damage to the limit of 
the variation among samples. 

On the other hand, if the recovery at -5°C is much weaker 
or slower than 25°C, the recovery during the four days of ir- 
radiation will reduce the damage to the 25°C group more. The 
results of the annealing at different temperatures indeed confirm 
this hypothesis and the recovery at higher temperature is faster 
and stronger. 



Figure 10: (Color online) The recovery of the signal amplitude. Signals from 
both samples returned to their levels before neutron irradiation. 
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Figure 1 1 : (Color online) Steps of the irradiation/annealing temperature depen- 
dence test. Numbers in boxes indicate number of samples tested in each group. 
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2. Right after the first irradiation, all the units were immedi- 
ately stored at three different temperatures, — 5°C, 25°C and 
60°C, for their first annealing. Every group had two units 
from each temperature group of the first irradiation. With all 
the units annealed, half of the units from the -5°C and 25 °C 
annealing groups were further heated to 60°C to investigate 
any additional recovery while the rest of the units were still 
kept at their original temperatures. 

3. After the first annealing was completed, all the units were 
irradiated again at -5°C or 25°C for four more days with an 
additional fiuence of 1.7 x 10 9 n eq /cm 2 , to see whether the 
radiation damage due to the first irradiation would effect the 
subsequent damage rate. 

4. At the end, all the units were heated to 40°C for a final ac- 
celerated annealing. 

4.2. Results 

4.2.1. Temperature Dependence of Radiation Damage and Re- 
covery 

The average current of all units before the irradiation test is 
86 + 3 nA and the uncertainty is the standard deviation of the 
measurements of individual units. After the first irradiation, the 
average current of the group at -5°C increased to 771 + 66 nA, 
and for the 25°C group, the current went up to 660 + 38 nA. 



Figure 12: (Color online) The recovery of the dark current at different temper- 
atures after the first irradiation. All the values were measured at 25°C and the 
uncertainties only include the accuracy of the measurements. 
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Figure 13: The secondary recovery at 60° C at the end of the first annealing of 
selected units previously annealed at -5°C or 25°C. 
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Figure 14: The final recovery of SiPMs at 40°C after the second irradiation. 



The recovery curves at various temperatures are shown in 



Fig. 12 13 and 14 The uncertainties in these plots only include 
the estimated uncertainty of the current measurement, 15 nA, 
while the variation among individual units is not included since 
it has no impact on the fit of the recovery time constant. The 
data are fitted with a simple exponential decay function with a 
constant baseline offset: 



I = b + a 



-t/T 



(1) 



where b is the baseline dark current after annealing, a is the 
recoverable damage and r is the time constant. 



Fig. 12 shows the recovery of the SiPM dark current at -5°C, 
25°C and 60°C right after the first irradiation. It is very clear 
that the units annealed at higher temperature recover faster and 
reach a lower asymptotic dark current. In order to see whether 
the baseline will be fixed after first annealing at a lower temper- 
ature, half of the units from the groups of -5°C and 25°C were 



heated to 60°C for a secondary annealing. As shown in Fig. 13 
additional recovery was observed and the time constant is con- 
sistent with the one obtained from the original 60°C group. The 
time constant of 40°C annealing were measured later after the 
second irradiation, as shown in Fig. [14] 

The temperature dependence of the recovery time constants 
can be described by an exponential curve 



t(T) = 41 • e 



-° lor day 



(2) 2 



is as plotted in Fig. 15 



Fig. [T6] shows the temperature dependence of the baseline. 236 
The uncertainties in the plot now include the variations from237 
individual units. Given the limited accuracy and the numbers 
of data points, the function of the temperature dependence can239 
not be well determined. It is clear nevertheless that when the24o 
annealing temperature is above 40° C, the change of the baseline24i 
can not be clearly identified given the variation of individual242 
units. 243 

As already discussed, the damage will further recover if244 
higher temperature is applied later. In order to see whether the245 



Figure 15: The dependence of the time constant T in Eq. l[TJ of the annealing 
on the temperature. Such a dependence is fitted by an exponential curve. 
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Figure 16: The dependence of the baseline b in Eq. Q of the annealing on the 
temperature. 



recovery would reverse when the annealed units are stored at 
a lower temperature, half of the units annealed at 60°C were 
put into a freezer for several weeks, and no indication of any 
increase in the dark current was found. 

The two units which were always kept at -5°C during the 
annealing after the first irradiation were annealed at 40°C with 
the rest of the units after the second irradiation. At the end, both 
units recovered to a level consistent with all the other units. 
This fact suggests that the temporary damage resulting from 
previous irradiations can always be recovered with sufficiently 
high annealing temperature. 

Finally, we plot the average dark current with annealing at 
temperature above 40° C as a function of neutron fluence in 



Fig. 17 The error bars represent the variation of dark current 
among units. The slopes of the damage during the two irradi- 
ations are consistent and it is clear that the previous irradiation 
will not effect later ones. 
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Figure 17: Damage curve of 1 mm SiPM as a function of 1 MeV neutron fiuence 
assuming annealing at 40-60°C. The current was measured at 25°C with a gain 
of'0.84xl0 6 . 



Figure 18: A typical ADC spectrum for the dark rate measurement. The data 
were taken from SiPM #1854 at 25°C after the first irradiation. The description 
of the fitting function can be found in the text. 



4.2.2. Relation between Dark Rate and Dark Current and Their 
Temperature Dependence 
What was being measured all the time is the dark current, 
but it is the dark rate which actually affects the performance of 
the Barrel Calorimeter in Hall D. In order to measure the dark 
rate, a DAQ system using a gated ADC was set up to record the 
random dark pulses generated by SiPMs. The length of the gate 
is chosen to be 200 ns which is longer than the pulse width of267 
SiPMs, -80 ns from 10% to 10%. Three 1 x 1 mm 2 SiPMs 0^ 
were measured at various temperatures between -5°C and 25°C269 
while the gain was kept constant by adjusting the bias voltagez™ 
based on the characteristic curve provided by Hamamatsu ' ifTTI .271 
In other words, the voltage setting above breakdown or over-272 
bias was kept constant. All three units were also part of thezra 
irradiation test therefore their ADC spectra after the irradiation™ 
were also taken for comparison. 275 
shows a typical ADC spectrum from a 1 x 1 mm 2 276 



Fig. 
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SiPM. The histogram was fitted by a convolution of a discrete277 
distribution function and a gaussian function. 27s 

279 

A(x) = Ao • P(n\p, Ap) ® Gaus(x\n ■ G + ped, cr„) (3) 280 

262 where Ao is the normalization factor. The discrete distribution, 281 

263 P(n\p, Ap), represents the probability that the number of pixels 282 

264 fired is equal to n, and it contains two Poisson distributions with 283 

265 one for the primary pixels fired and the other for the total of 

266 cross talk or after pulses 1 12 1 caused by the primary pixel: 

P(n\p, Ap) = Pois(i\p) ■ Pois(j\i ■ Ap) 



n=i+j 



,-(/i+/A/i) / 



y{iApy 



(4) 



where p is the average number of primary pixels fired and Ap is 2 
the average number of pixels fired around the primary pixel due 



to cross talk and after pulses. The Gaussian function, Gaus(x\n ■ 
G + ped, cr„),represents the distribution of the charge with the 
number of pixels fired equal to n. G is the total gain |^] in ADC 
channels and ped is the ADC pedestal value. The width, <x„, is 
equal to 

o-„ - Jt 



jcr 2 ,+n-cr 2 

ped sig 



(5) 



where <r pe d is the width of the pedestal and cr J(g is the intrinsic 
width of a single pixel signal. Such a function is valid when 
both the signal occupancy and the Ap are small, which is true 
for our test condition. Otherwise, the function needs to be mod- 
ified by replacing the Poisson distributions with binomial dis- 
tributions. 

The absolute gain of SiPMs was calculated by dividing the 
total gain G by the known ADC conversion factor and the gain 



of the preamplifiers. As shown in Fig. 19 the gain is relatively 
stable at different temperatures as long as the bias voltage set- 
ting was adjusted to compensate for the change in the break 
down voltage with temperature to have a fixed over-bias. The 
uncertainties shown in the plot only include the statistical un- 
certainty from the fit and the uncertainty of the temperature 
reading. However, due to the temperature gradient in the cool- 
ing device, the gain fluctuates slightly around the average value, 
0.84 x 10 6 . 

Fig. [20] shows the correlation between the dark current and 
the dark rate, and the dark current has been corrected for the 
deviation of the gain to the average value. The correlation is 
well described by a linear function: 



/ = 0.190 nA/kHz ■ / + 3.84 nA 



(6) 



Their serial numbers are 1853, 1854 and 1855. 



The slope corresponds to an average gain of 1.19 x 10 6 which 
is about 42% higher than the actual gain. Part of the mis- 
match comes from the cross talk and after pulses which are 

6 It includes the intrinsic gain of an APD pixel, the gain of the pre-amplifiers 
and the ADC's analog-to-digital conversion factor. 
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Figure 19: The absolute gain of SiPMs at different temperatures at constant 
over-bias. Different legends represent the data from different units before and 
after the irradiation. 



Figure 21: The dependence of cross talk on temperature. 
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Figure 20: The correlation between the dark current and dark rate with a gain 
of 0.84 X 10 6 . 



not counted in the extraction of the dark rate. The rest may 
be attributed to the fact that Eqn. <[3j does not fully account for 
the pulses partially integrated in the ADC gate Q Clearly, the 1 
radiation damage does not have any impact on this relation. 

On the other hand, there is no significant temperature de- 
pendence of the cross talk and after pulses observed as shown 
in Fig 
well. 



21 And the radiation damage doesn't change them as 



Fig. 22 shows the temperature dependence of the dark rate?^ 
before and after the neutron irradiation. The dependence is ex-3°° 
ponential in the measured temperature range and the change of™ 
the temperature coefficient caused by the radiation damage isso2 
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7 In retrospect, we realize that a longer gate width may have allowed a more 1 
accurate determination of the extracted dark rate from the fit. 
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Figure 22: The dependence of the dark rate before and after irradiation on 
temperature. 



relatively small. As a result, the average dependence is 

/ = /o . e 0.075.(T-7-„) 



(7) 



Such a behavior provides the motivation to cool SiPMs during 
beam time in Hall D to reduce the dark rate. 



297 5. Summary 



We measured the neutron radiation damage to SiPMs using 
neutrons generated by an electron beam at Jefferson Lab and a 
AmBe neutron source. We further studied the temperature de- 
pendence of the radiation damage and other properties includ- 
ing dark rate, dark current and damage recovery. We found that 
both dark rate and dark current increase linearly as a function 
of the total neutron fiuence and the damage does not depend on 
the temperature or operating voltage. Part of the acute damage 
will recover. The speed and the extent of this annealing process 



8 



307 strongly depends on the temperature and is faster and stronger 

308 at higher temperature. Increasing the temperature of a dam- 

309 aged unit previously annealed at a certain temperature brings 

310 further recovery, but lowering the temperature will not reverse 

311 the recovery achieved. We also measured the temperature de- 

312 pendence of the dark current and dark rate of SiPMs at a fixed 

313 gain. Such a dependence is not strongly affected by the neutron 

314 radiation damage. 

315 The results obtained by this study provided important infor- 

316 mation for implementing SiPMs as the readout of the Barrel 

317 Calorimeter in JLab Hall D. 
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